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ABSTRACT 

The first population III stars are predicted to form in minihalos at z w 10-30. 
The James Webb Space Telescope (JWST), tentatively scheduled for launch in 2018, 
will probably be able to detect some of the first galaxies, but whether it will also be 
able to detect the first stars remains more doubtful. Here, we explore the prospects of 
detecting an isolated population III star or a small cluster of population III stars down 
to z = 2 in cither lensed or unlensed fields. Our calculations are based on realistic 
stellar atmospheres and take into account the potential flux contribution from the 
surrounding H n region. We find that unlensed population III stars are beyond the 
reach of JWST, and that even lensed population III stars will be extremely difficult 
to detect. However, the main problem with the latter approach is not necessarily 
that the lensed stars are too faint, but that their surface number densities are too 
low. To detect even one 60 M population III star when pointing JWST through 
the galaxy cluster MACS J0717. 5+3745, the lensing cluster with the largest Einstein 
radius detected so far, the cosmic star formation rate of population III stars would 
need to be approximately an order of magnitude higher than predicted by the most 
optimistic current models. 
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1 INTRODUCTION 

Both theoretical argum ents and numerical simulations (e.g. 
iBromm fc Larsonl [2004 ) , strongly support the notion that 
population III stars were very massive, significantly more so 
than the population I and population II stars that formed 
later on. When stars form in metal-enriched gas, the Jeans 
mass is lower, which leads to fragmentation and thus lower 
masses. Lacking metals, a chemically unenriched cloud does 
not fragment in the same way and the characteristic stel- 
lar mass may therefore be higher. It has been argued that 
two different classes may have existed: population III. 1 and 
population III. 2. Population III. 1 stars formed in dark mat- 
ter minihalos of mass 10 5 -10 6 M Q at z g 10-30. Early re- 
searc h on population HI star formation dO'Shea, fc NormarJ 
120071 ; lYoshida et ail 120081 : IBromm et al.l l2009r i. indicated 
that only one star, with a very high average mass of ~ 100 
Mq was produced in every minihalo. This was the result 
of UV radiation (in the Lyman- Werner band) produced by 
the first massive star destroying the molecular hydrogen in 
the parent cloud, preventing further cooling and star forma- 



tion. Later research on the same t opic (e.glKrumholz et ahl 
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Greif et al.ll2011l : 

cates a more complex scenario with more fragmentation, re- 
sulting in a binary system or small clust er of stars, with 
the m ost massive star reaching f» 50 Mq . IStacv fc Bromnj 
(2008) also suggest that cosmic rays from the supernovae re- 
sulting from the first massive population III stars could sig- 
nificantly influence subsequent star formation in minihalos, 
possibly yielding a characteristic stellar mass of ~ 10 Mq. 
It is also plausible that protostellar feedback will halt the 
mass accretio n, thereby limiting th e stellar mass to a few 
times 10 M (|Hosokawa et aUbOllT ). 

There are also models indicating late, 2 < z < 7, pop- 
ulatio n III star formation jTornatore et al.l 120071 : Ijohnsor] 
2010). These late population III stars form in regions 
that are pristine due to highly inhomohegeneous metal 
enrichment. These low-metallicity regions are usually lo- 
cated in the vicinity of more massive (~ 10 s — 10 11 M©) 
dark matte r halos containing metal enriched galaxies. In 
llnoue et al.l (|201ll 'l. the authors examine LAEs at z = 
3.1 with signatures they argue are hard to interpret as 
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other than objects contain i ng po pulation III stars, (but see 
iForero- Romero &: Diikstral (|2012h for a different view). 

Population III. 2 stars are thought to form mostly due to 
cooling by adiabatic expansion or HD cooling, usually at a 
lower redshift than population III. 1 stars, and are thought to 
be building blocks for early galaxies. Since we focus on pre- 
galactic population III star formation we will not investigate 
thes e any further. 

iBromm et al.l (|200ll ) have investi gated the spectral en - 
ergy distribution of a p rimordial star. iGardner et ail (120061 ) 
and iGreif et al.l (|2009l ) used these results to show that iso- 
lated population III stars are too faint for detection with the 
JWST. However, the flux boost due to gravitational lensing 
has not previously been considered in detail, e xcept for a pre- 
limin ary version of the study we present here (|Rvdberg et al.l 
l2010l ). 

Apart from direct detection through gravitational lens- 
ing, there are other ways to probe the properties of popu- 
lation III stars. The cumulative indirect signatures in the 
global radio background produced by bremsstrahlung or 
the 21-cm hydrogen emission also represent interesting op- 
tions. Bremsstrahlung is free-free emission originating in the 
Hn region in its active and relic states. When the popula- 
tion III star has died the almost fully ionized Hn region 
around it begins to recombine. This makes it a very brigh t 
source of hydroge n 21-cm radiation (jTokutani et alj 12009). 
IGreif et al. (2009) find that bremsstrahlung is very unlikely 
to be detected in the near future but the cumulative 21- 
cm radiation might be strong enough to be within reach of 
the upcoming Square Kilometer Array. Direct observation 
of pair instability supernovae (PISN) could also be used to 
infer the existence of population III stars. PISNs are the 
predicted fin al fate of non-rotat i ng stars in the mass range 
140-260 Mm dBarkat et al.ll 19671 ; lFralevlll968l ; IWhalen et all 
120121 ; iDessart et al.ll2012h. and are possible for lower mass 
stars if rotating l|Chatzopoulos fc W heeler 2 0121 ). The initial 
mass function (IMF) describes the distribution of mass for 
a new star. Since PISN occur in a certain mass range obser- 
vations of those could probe the IMF. Using a sem i-analytic 
halo mass function approach. iHummel et al.l (|2012| ) find that 
no more than ~ 0.2 population III PISN should be visible 
per JWST field of view at any one time, depending on the 
adopted model for stellar feedback. Rapidly rotating popula- 
tion III st ars may also produce collapsar gamma-ray bursts, 
or GRBs l|de Souza et al.ll201ll ). These are formed by metal 
poor, rotating stars more massive than ~ 40 Mq . This could 
be a highly energetic probe of population III stars, and of 
the high-mass tail of their IM F. Population III stars could 
also leave detectable imprints (|Kashlinskvll200fJ) in the cos- 
mic infrared background (CIB). The CIB is a repository of 
emissions from different sources at different redshifts, all red- 
shifted into the infrared. The imprints from population III 
stars are already cons trained by det ection of emission at en- 
ergies above 10 GeV (|Gilmorell2012T ). as interaction between 
the radiation of population III stars and gamma rays create 
a gamma-ray optical depth. 

As we will argue, a flux boost due to gravitational lens- 
ing by a foreground galaxy cluster will be necessary to di- 
rectly detect population III stars. By adopting the properties 
of one of the more promising lensing clusters, we will calcu- 
late the SFR required to achieve ~ 1 detectable population 
III star per survey field. 



2 MODELING THE SPECTRA OF 
POPULATION III STARS 

We have used the p roperties of population III stars derived 
by ISchaereil l|2002l ) as the basis for our estimated broad- 
band fluxes. Here, we adopt the main sequence properties of 
these stars and thereby neglect the fact that their ionizing 
fluxes decrease as they grow older. This results in a mild 
overestimate of the flux in the filters we are using. We have 
used realistic stellar atmospheres and have also taken into 
account the potential flux contribution from the surround- 
ing Hn region. To model the stellar atmo spheres, we have 
used the publicly available TLUSTY code, |Hubenv fc Land 
1995). This code computes ID, non-LTE, plane-parallel stel- 
lar model atmospheres and spectra, given a certain input 
stellar composition, surface temperature and gravity. We set 
the chemical comp osition of our star s to the primordial mix- 
ture of H and He (|locco et al.ll2009l ). 

For our models of population III stars, the luminosity 
and number of ionizing photons (important for the lumi- 
nosity originating in the nebula around the stars) scale ap- 
proximately linearly with stellar mass for M > 25 Mq, but 
not below this limit. This means that the luminosity will be 
dominated by a few mass ive stars, one or possibly two of a 
few tens of solar masse s (Krumholz ct al. 2009; Tur k et ail 
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Their combined luminosity will correspond to a star with 
stellar mass slightly lower than the sum of their masses, due 
to the non-linearity of the luminosity as a function of stellar 
mass below 25 Mq. Smaller stars add very little compared 
to larger stars. For example, a 9 Mq star has an energy out- 
put (bolometric and ionizing) that corresponds to less than 
10% of that of a 25 Mq star. In this paper we will work with 
one population III star as a proxy for small clusters of this 
type. We select a 60 Mq star as our default option because 
its output likely exceeds the output of a smaller cluster. As 
an example of a more extreme scenario, we will also present 
calculations for a stellar mass of 300 Mq; this is the lowest 
stellar mass that turns out to be readily detectable in our 
models. A 200 Mq star, for example, is only detectable for 
a small subset of models, and only over a small redshift in- 
terval. Our proxy stars admittedly have lower lifetimes than 
the constituents of the cluster they are representing. This 
means an argument could be made for using longer lifetime 
in our calculation, thereby improving the prospects of detec- 
tion. Because the lifetimes do not differ by much (a 40 Mq 
and 60 Mq star differ 10% in lifetimes, as an example), and 
it would be difficult to make a more accurate estimate of 
the effect, we will adopt the lifetime of the proxy star as 
representative of the lifetime of a small cluster. 

The observed flux from a population III star is prob- 
ably st rongly influenced by the Hn region surround- 
ing it (|Kitavama et al.l l2004h . To model this, we have 
used the publicly a vailable photoionization code Cloudy 
l|Ferland et al.| [l998). Assuming spherically symmetric Hn 
regions, we use Cloudy to calculate the resulting spectrum. 
We predict both the nebular continuum and emission lines. 
This procedure results in an optimistically bright emis- 
sion spectrum, as more realistic nebulae experience feed- 
back effects that could give rise to holes in the Hn re- 
gion. This exposes the star underneath and the spectrum 



Detection of population III stars with the JWST 3 



emerging from the stellar atmosphere. The region could 
also break out of the gravitational well of the host halo if 
the feedback is stro ng enough (e.g. iKitavama et al.l 
IWhalen et"afl |2004 ). As described by iGreif et all 
this dilutes the nebular flux. This could expose more of the 
purely stellar continuum resu lting in a different overall spec- 
trum l|Zackrisson et al]|201ll ). 

For the low densities expected in the H n regions around 
the stars we have considered, time dependent effects start 
to become important at late times when the density may 
become as low as f atom cm -3 , or lower. As Cloudy is a 
steady-state code, we have tested the importance of time- 
dependence using the latest version of the time-dependent 
photo io nization code de s cribed in iLundqvist &: Franssonl 
\ 19961 ). iLundavisd l| 19991 . 120071 ) and iMattila et all l|20ld ) 
For hydrogen and helium, a total of 56 levels for H I— II and 
72 levels for Hel-lll are included, and all levels are treated 
time-dependently. We discuss the importance of time depen- 
dence for the emission from the Hn regio n in Section [H 

A t z > 6 the Gunn-Peterson trough l|Gunn fc Peterson! 
1 19651 ) affects each emitted spectrum, resulting in huge ab- 
sorption of radiation with wavelength lower than the Lyman- 
a (Lya) line. To simulate this we set all flux at wavelengths 
shorter than Lya to zero at z > 6. For z < 6 the situation 
is more complicated as the radiation passes through several 
clouds of neutral hydrogen. Each cloud absorb Lya radia- 
tion producing an absorption line. The resulting absorption 
lines are called the Lya-forest. The exact nature of the Lya- 
forest is dependent on the line of sight with the ave rage ab- 
sorpti on increasing with redshift. We have used the lMadarJ 
1 19951 ) model to modify the radiation below the Lya-line 
at z < 6. Even though Lya is the main absorbing mech- 
anism the model also takes higher order Lyman-absorbers 
into account. We use the average absorption predicted by 
the model. 

The very high effective temperature of a population III 
star (~ 100,000 K) implies a huge hydrogen ionizing flux. 
This is converted into lower energy radiation but also into 
a substantial number of Lya photons. This can potentially 
make the Lya line the strongest rest frame emission line in 
the spectra of these population III stars. At the same time, 
this is a highly resonant line which means that it suffers from 
absorption in the neutral intergalactic medium (IGM). The 
Lya photons must also escape the interstellar medium (ISM) 
where they scatter and might be destroyed by dust. Even so, 
if just a fraction of these photons manage to escape the ISM 
and IGM absorption they could significantly improve the 
prospects of detecting population III stars. The combined 
fraction of Lya photons escaping both the ISM and the IGM 
will be denoted fhya- Our default assumption is to set the 
Lya line to zero (/Lya = 0), but we will also consider two 
scenarios with strong escape fraction, setting fhya to 0.2 or 
0.5. 



60 M, 




Figure 1. The NIRCam/F200W AB magnitude predicted for a 
60 Mq star as a function of redshift in an unlensed field. The blue 
line represents the magnitude with no Lya radiation escaping, 
whereas the purple and black lines represent /Lya = 0.2 and 0.5, 
respectively. The red dashed line indicates the detection limit 
(30.9 mag) considering S/N = 10, using an exposure time of 
lOOh. As seen, the unlensed flux of a 60 Mq population III star 
lies > 6 magnitudes below the JWST detection threshold. 



detection of population III stars at the redshifts we are in- 
terested in. The brightest population III magnitude is 36.6 
at a redshift of 2. For comparison, the S/N = 10, lOOh ex- 
posure time detection limit is 30.9, placing this population 
III star > 6 magnitudes below the JWST detection limit. 

As it appears impossible to detect a high redshift pop- 
ulation III star in an unlensed field, we will in the follow- 
ing instead explore the prospects of detecting population 
III stars in gravitationally magnified regions of the sky. For 
this purpose we adopt MACS J0717. 5+3745, a galaxy clus- 
ter at z — 0.546, as our lensing cluster. This is the gravita- 
tional lens with largest angular Einstein radius known, and 
also one of the t argets of the ongoing HST survey CLASlR 
jPostman ct al. 2012|). 

IZitrin et al.l (|2009l ) presented a detailed model of this 
particular gravitational lens. The unusual structure of the 
object, with an unrelaxed morphology and correspondingly 
shallow density profile, results in exceptionally large areas 
with high magnification. For instance, at z > 6, the area over 
which background objects achieve a magnification higher 
than 10 is w 3.5 arcmin 2 . The data we have for the mag- 
nification of MACS J0717.5+3745 cover z > 6. For lower 
redshifts we have used the data for z — 6, this should only 
have a small impact on our calculations as the magnified 
area changes slowly with redshift. 



2.1 MACS J0717.5+3745 

To assess the prospects for directly detecting a 60 Mq pop- 
ulation III star, we calculate the intrinsic JWST broadband 
flux of such a star as a function of redshift and compare it 
to the detection limit. In Figure [T] we display the results for 
the JWST/NIRCam F200W filter. The sensitivity of NIR- 
Cam when using F200Wis generally the most promising for 



3 THE SFR REQUIRED FOR JWST 

DETECTION OF POPULATION III STARS 
IN LENSED FIELDS 

Figure [2] shows the minimum magnification (/i m i n ) required 
for detection of one 60 Mq star as a function of redshift in 
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the J0717.5+3745 field. At each redshift the filter requir- 
ing the lowest magnification is used when calculating fi m - m . 



60 Mo 



We consider a generous exposure time of t e 



100 hours 



and use a detection limit of S/N — 10 per pointing^. For 
a fixed detection limit, S/N oc ^/£ e xp- This means, as an 
example, that all results hold for S/N = 2 using i cxp = 4h 
and for S/N = 5 using t oxp = 25h, as well as our default 
S/N = 10 using the t cxp = lOOh case. Included are also two 
different non-zero Lya escape fractions: ft,ya = 0.2 and 0.5. 
A positive fh ya facilitates detection at most redshifts. For 
redshifts where the Lya line is located between filters the 
required magnification is of course unchanged in each filter. 
For a 60 Mq star, magnifications > 200 are required, as seen 
in Figure [2] Magnifications this large are probably not im- 
possible, but will b e attained only ov er very small areas in 
the lensing cluster (|Zitrin et al.ll2009l ). 

Using the result from Figure [2] we calculate the to- 
tal lensed area behind MACS J0717.5+3745 that displays 
the required magnification for detection in a unit red- 
shift interval as a function of redshift. The correspond- 
ing co-moving cosmological volume is calculated using 
Ho=72 km s" 1 Mpc -1 , Q m = 0.27 and Q A = 0.73. We 
have adopted t he lifetimes of population III stars given in 
ISchaererl |2002) for this calculation, and considered the star 
formation rate (SFR) constant throughout this time interval 
at each redshift. By combining the adopted typical mass of 
population III stars with their predicted lifetimes, we arrive 
at a lower limit (SFR m i n ) on the SFR required to form at 
least one population III star within this volume. 

We have explored thre e differ ent star formation mod- 
els from ITrenti fc Stiavellil i2009h. the simulations from 
IWise et all (|2012h and iTornatore et all (I2007T). to compare 
with our results for SFR m i n . The three ITrenti fc Stiavellil 
l|2009l ) models use Eb cooling in minihalos to form the stars, 
the difference between them being how they handle the ex- 
pected Lyman- Werner (L-W) radiation feedback. L-W radi- 
ation is emitted by hot stars and destroys H2 through pho- 
todissociation, thereby prohibiting further population III 
star formation. The first, standard, model assumes full feed- 
back from L-W which effectively prevents population III star 
formation at z < 13. The other two models assume reduced 
and no L-W feedback. These scenarios correspond to situa- 
tions where the L-W radiation is lowered for some reason. 
This could happen in the case where the "isolated" popula- 
tion III stars we consider in fact belong to small clusters with 
many low-mass stars with low L-W fluxes, thereby decreas- 
ing the effective L-W flux per unit stellar mass. Instabilities 
in the ionization fronts resulting in self shielding could also 
ha mper L-W fe e dback (IWhalen fc Normanl [20081 ). 

IWise et all (|2012l ) use an adaptive mesh radiation hy- 
drodynamics simulation to model star formation. In their 
models, the ionizing and dissociating radiation from nearby 
population III stars and galaxies proves to be unable to 
halt population III star formation. Instead, it only sup- 
presses it, resulting in a nearly constant SFR around 5 x 
10~ 5 Mq yr -1 Mpc~ 3 throughout the redshift interval z = 
7-18. 



Note that it is possible that it takes two pointings to 
cover the J0717. 5+3745 caustic (high magnification area) with 
JWST/NIRCam. 




Figure 2. The minimum gravitational magnification /i m in re- 
quired to detect a 60 Mq population III star as a function of 
redshift. The detection criterion is based on a 10cr detection after 
an exposure time of lOOh. For each redshift, the filter resulting in 
the lowest necessary magnification is selected in the calculation 
of Mmin- Clearly, direct detection without gravitational lensing is 
impossible and even when seen through a lensing cluster, a very 
high magnification (/^ > 200, increasing steeply with redshift) is 
required. The blue line represents the lowest required fi with no 
Lya radiation escaping, whereas the purple and black lines rep- 
resent /Lya = 0.2 and 0.5, respectively. 



To compare our results to lower redshift population 
III star formation, possible through inhom ogeneous metal 
enrich ment, we use the standard model in ITornatore et all 
|2007T ). Their code uses smoothed particle hydrodynamics to 
simulate the development of star formation. The suppression 
of star formation at z < 6 is caused by the IGM photoheat- 
ing due to reionization. The collapse of low-mass structures 
is inhibited as the photoheating increases the Jeans scale. 

In Figure we plot SFR m j n against redshift for the 
case of a 60 Mq star. This graph includes the two cases 
with non-zero fhya and our comparison SFR models. There 
is a gap of an order of a magnitude between the most op- 
timistic simulations and SFR m i n , even when considering a 
very high Lya escape fraction. Thus, our default scenario is 
undetectable. The corresponding requirements for a 300 Mq 
star are shown in Figure |U In this detection is possi- 

ble if considering one of the most optimistic SFR predictions 
and a very high Lya escape fraction. 

The IMF of population III stars is essentially unknown. 
Simulations indicate a high typical mass but the exact dis- 
tribution is important for the prospects of detection. In Fig- 
ure SFR m i n is plotted as a function of the typical mass 
of population III stars at z — 15 in the filter F200W. We 
selected this redshift and filter passband in order to mini- 
mize SFRmin, while simultaneously maximizing the poten- 
tial impact of the Lya emission line. To produce a given 
number of stars, higher stellar masses imply a higher SFR, 
in order to supply the additional mass. More massive stars 
also live shorter lives, further increasing SFR m in- On the 
other hand, a higher mass implies a higher luminosity. This 
means that a larger cosmological volume in a unit redshift 
has sufficient magnification for observation, implying a lower 
SFRmin- It turns out that the effect of the luminosity trumps 
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Figure 3. The minimum star formation rate SFR m i n required 
for approximately one 60 Mq star to be detected when JWST 
is pointed through the galaxy cluster MACS J0717. 5+3745. The 
detection criterion is based on a 10c detection after an exposure 



time of lOOh. The blue solid line represents the SFR n 



fith 



no Lya radiation escaping, whereas the purple and black solid 
lines represent /Lya = 0.2 and 0.5, respectively. Different SFR 
models are also included for comp a rison, using dashed lines. The 
SFR models by iTrenti fe Stiavellil ||2009|) are represented by the 
green, orange, and red lines. They represent different levels of L-W 
radiation feedback, a feedback prohibiting further star formation. 
The three levels are basic (green), meaning full L- W feedback, 
reduce d (o range), and no L - W (re d). The models by I Wise et al.l 
||2012| 1 and iTornatore et al. 1 11200711 are included in cyan and grey 
respectively. The figure shows that there is a gap of an order 
of magnitude between the most optimistic simulation and the 
required SFR for a detection, even when a very high Lya escape 
fraction is assumed. 
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Figure 4. Same as Figure [3] but for a 300 Mq star. With this 
fairly extreme stell ar mass, a detection wou ld be possible with 
the most optimis tic ITrenti fe Stiavellil ||2009| ) simulation and the 
I Wise et all {20121 ) model, but only when a very high Lya escape 
fraction is assumed. The mass 300 Mq is selected for the star 
to be readily detectable with the star formation models we com- 
pare to. A 200 Mb star , for example, is not detectable with the 
ITrenti fc S tiavclli ( 200|) models even though it st i ll is d etectable 
at some redshifts when employing the lWise et al. I HoTJ) model. 



Figure 5. The SFR min in the MACS J0717.5+3745 field, as a 
function of the population III mass (assuming all such stars to 
have the same mass) at 2 = 15. Imaging in the NIRCam/F200W 
filter is assumed. Since this graph uses a single filter, it is not 
necessarily the filter resulting in the lowest SFR m ; n . This means 
that the values in this graph differ in some instances from the 
values in Figure [3] and [4] where the filter resulting in the lowest 
SFR m ; n is always selected. The detection criterion is based on a 
lOcr detection after an exposure time of lOOh. The blue line rep- 
resents the SFR m ; n with no Lya radiation escaping, whereas the 
purple and black lines represent /Lya = 0.2 and 0.5, respectively. 
A higher population III stellar mass results in a lower required 
SFR. This indicates that a top heavy IMF improves the prospects 
of detecting high-redshift population III stars. 



the mass and lifetime effects. Hence, a larger mass implies a 
lower required SFR, as can be seen in Figure [5] Thus a top- 
heavy IMF improves the prospects of detecting population 
III stars. 

So far, we have concentrated on the prospect of detect- 
ing population III stars in a specific redshift bin by selecting 
the filter resulting in the lowest SFR m i n . Broadband imaging 
surveys will in reality probe a wide range of redshifts simul- 
taneously. For the purposes of comparison, here we consider 
the redshift interv al 7-15, which roughly corresponds to the 
IWise et al.l (|20 12f ) model. For computational purposes, we 
kept the SFR constant in the redshift interval. We arrive at 
a slightly lower cosmic SFR m i n compared to the redshift bin 
appro ach. This is more or less consistent with the lWise et alj 
(2012) model, but in reality could of course result in a more 
complicated SFR. Table[JJdisplays the result for the different 
JWST filters in the NIRCam and MIRI instruments. Com- 
paring to the lowest values in Figure [3] we see that larger 
redshift windows do not significantly alter our conclusions. 
This conclusion follows mostly from the IGM absorption. In 
the filters that return the lowest SFR m i n (which occurs at 
low redshifts), the radiation coming from higher redshifts is 
mostly absorbed by the IGM. Hence, the final SFR m i n in 
those filters is dominated by a small interval at low redshift. 



4 DISCUSSION 

Of the many scenarios explored in Section [31 the only one 
that would make isolated population III stars potentially 
detectable with JWST is: 
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Table 1. The minimum logio SFR m ; n (Mq Mpc — 3 yr — 1 ) required for detection of one 60 Mq population 
III star anywhere within the redshift interval 7 < z < 15, as a function of the filter used. A constant SFR 
over the relevant redshift interval is assumed. The detection criterion is based on a 10o" detection after an 
exposure time of 100 h. Three cases are considered: no escape through ISM and IGM of Lya photons and 
the cases /Lya = 0.2 and 0.5. The best filter choices for each /Lya have been marked in boldface. The 
radiation through the lower wavelength filters F070W and F090W are to the most part absorbed by the 
Gunn-Pctcrson trough and so are excluded. We do not include MIRI for the longer wavelength filters, as its 
sensitivity is too low. 



Filter 


Instrument 


A (pm) 


t Lya = 0.0 


tjjya = 0.2 


fiya = 0.5 


F115W 


N IRC am 


1.15 


-2.38 


-2.95 


-3.44 


F150W 


N IRC am 


1.50 


-2.99 


-3.14 


-3.41 


F200W 


NIRCam 


1.99 


-3.23 


-3.27 


-3.36 


F277W 


NIRCam 


2.79 


-3.04 


-3.04 


-3.04 


F356W 


NIRCam 


3.56 


-3.03 


-3.03 


-3.03 


F444W 


NIRCam 


4.44 


-2.64 


-2.64 


-2.64 


F560W 


MIRI 


5.63 


-1.19 


-1.19 


-1.19 


F770W 


MIRI 


7.65 


-0.72 


-0.72 


-0.72 



• Very massive (> 300 Mq) population III stars. 

• The stars are formed in the limit ing case of either no 
L-W feedback (|Trenti fc Stiavellill2009h or in th e prolonged 
popul ation III star formation case predicted bv lWise et al.l 
|2012l ). 

• The Lyman-alpha escape fraction would need to be very 

high (/ L ya w 0.5). 

• The gravitational lens contains sufficiently large regions 
with very high magnifications (> 1000). 

• The escape of Lyman continuum (i.e. hydrogen- 
ionizing) photons from the nebula surrounding the popu- 
lation III stars is very small. 

Hence, the prospects for detecting isolated population 
III stars in the foreseeable future appear bleak. However, 
contr ary to previous claims (|Gardner et al.l l2006 ; Grci f et all 
2009), the problem is not necessarily that population III 
stars are too faint for detection, as the magnification of 
a suitably chosen lensing cluster may lift them above the 
JWST detection threshold. The main obstacle is instead that 
the surface number densities for sufficiently massive popu- 
lation III stars to be found in a sufficiently lensed field is 
likely to be too low. 

Even so, given the scenario that we have a detection 
of a population III star we still have to identify it as such. 
One plausible mean of doing this suggested in the litera- 
ture (|Oh et all l200ll ; iTumlinson et al.ll200lh is to identify 
the Hell 1640 A line. Population III stars are predicted to 
have a significantly harder spectrum than even low metal- 
licity stars. This means significantly more He are ionized 
in relation to H ionized. The result is that the ratio of the 
Hell 1640 A line to a non-resonant hydrogen line such as 
Ha should be significantly larger for population III stars 
compared to that from even low metallicity stars. 

It has also been argued that supersonic streaming ve- 
locities between dark matter and gas could delay popula - 
tion III star formation, (|Greif et alj|201ll ; lMaio et al-lbOllT ). 
The delay is substantial, Az ~ 4, but the large shift in 
mass also means the number density of population III star- 
forming minihalos is reduced by up to an order o f magnitude. 
Ther e is criticism of the importance of the effect (|Stacv et al.l 
l201ll ), but if significant, it would of course impact the re- 



sults presented here. The combined effect of fewer stars at 
lower redshift needs careful examination in future work. 

As mentioned in Section[2] the low densities expected in 
the Hn regions could make time dependence important, in 
particular when the time scales of ionization/recombination 
and cooling/heating become comparable to the lifetimes of 
the stars. We have tested this using two scenarios for the 
surrounding H n regions, one with 1 and another with 100 
hydrogen atoms cm -3 . We set the He/H number ratio to 
0.09, and the abundances of other elements to zero. We show 
the relative luminosity of Lya for the two densities and for 
both the 60 M Q and 300 M Q stars in Figure [6] We have 
assumed that the stars suddenly switch on and shine for 
about 3.4 Myrs with a constant luminosity and spectrum. 
For tih = 100 cm -3 , the ionization quickly goes toward the 
steady-state solution, whereas for nu — 1 cm" 3 it takes 
~ 0.4 Myrs to reach the steady-state level. However, even 
after this level has been reached, there is still some slow 
photoionization going on in the outer parts of the H II re- 
gion. The temperature in this region also overshoots some- 
what compared to the steady-state solution before steady 
state has been attained, which produces slightly enhanced 
collisional excitation of Lya. This effect is larger for the 
300 Mq star with its more extended region of partially ion- 
ized hydrogen outside the Stromgren radius. The net effect 
of time dependent versus steady-state models is, as can be 
seen from Figure [6l not more than ~ 10 — 15 %. For den- 
sities below riH = 1 cm" 3 , the deviations from steady-state 
will obviously become larger. However, the Hn region is ex- 
pect ed to be evolv i ng fro m high to low density with time 
(e.g. IWhalen et alj 12004'). This probably makes the early 
delay of Lya in Figure exaggerated as, at this time, the 
Hn region is dense. At low densities the emission from the 
H n region will fall off on a time scale comparable to the life- 
time of the star. For such low densities, the H n region may 
however no longer even be ionization bounded. The uncer- 
tainty in our models due to time dependence is in any case 
much smaller than due to other assumptions. 

We have calculated the required magnification for de- 
tection of population III stars using realistic models for the 
atmosphere and the Hn region. However, our results are ad- 
mittedly based on a number of simplifications, most of them 



Detection of population III stars with the JWST 7 



1 1 1 1 




■ 


60 M s 




■ 


300 M s 


■j L 


1 1 ■ ■ 


i 



1 2 

Time (Myr) 



Figure 6. The flux in the Lyman alpha line from a njr = 1 cm" 3 
nebula divided by the flux in the Lyman alpha line from a njj = 
100 cm" 3 nebula, as a function of time. We have assumed that 
the stars suddenly switch on and shine for about 3.4 Myrs with a 
constant luminosity and spectrum. For the reference nebula with 
nji = 100 cm -3 the ionization quickly goes toward the steady- 
state solution. The blue line represents a 60 Mq star, whereas the 
purple line represents a 300 Mq star. For njj = 1 cm" 3 it takes 
~ 0.4 Myrs to reach the steady-state level. A 60 Mq star's nebula 
subsequently emits ~ 5 % more flux than in steady state and a 
300 Mq star's nebula up to ~ 15 % more flux than in steady state. 
The upturn seen for both stellar models after 3.4 Myrs is an effect 
of time dependence in the low density (1 cm -3 ) case. For both 
densities, the Ly-alpha emission starts to decay when the star 
is switched off, but the decay is about two orders of magnitude 
faster for 100 cm" 3 than for 1 cm -3 . 



made in the direction that would improve the prospects of 
detection. The spherically symmetric Hn region that we 
have assumed likely overestimates the flux, as there proba- 
bly will be low-density channels in the nebula through which 
ionizing radiation would be able to escape into the inter- 
galactic medium. We have also assumed a spatially homoge- 
neous distribution of minihalos. In reality clustering could 
impact the detection probability in a negative way. All in all, 
even though not impossible, the prospect of detecting pop- 
ulation III stars (isolated, or in small clusters) with JWST 
appears bleak at best. 

Larger numbers of population III stars could poten- 
tially form in chemically p ristine (10 7 -10 8 Mq) atomic - 
cooling halos at z < 15 (e.g. Saf ranek-Shrader et al.ll2012h . 
These "population III galaxies" woul d by comparison be 
much easier to det e ct and also identify jjohnson et "al1l2009l ; 
iPawlik et all 1201 ll ; IZackrisson et all 120111 ). The detectabil- 
ity of such objects hinge more on the total star forming 
mass than on stellar masses and nebular properties, re- 
sulting in less ex t reme lensing magnification requirements 
|Zackrisson et all I2012T ). The merits of investigating the 
prospects for detecting population III stars in minihalos is 
that they possibly fill an important role in the development 
of the universe, providing the first metal-enric hment prior 
to the first galaxy formations l|Greif et al.ll201oT ). 

More speculative and exotic theories involving popula- 
tion III stars or derivatives of them also exist, potentially im- 
proving detection prospects. Fast accretion combined with 
rapid rotation could potentially form super massive popu- 



lation III stars (Hosokawa et al. 2012T). These would be ex- 
tremely massive, up to 10 5 -10 6 Mq, and have lower tem- 
perature, ~ 5000 K. A black hole could also form by col- 
lapse in the center of the star powering it through emis - 
sion by accretion, a so-called quasi-star |Begelmanll2010t ). 
Some significant fraction of population III stars could also 
be "dark stars" , which would make the detection prospects 
slightly better (Zackrisson et al. 2010a). There are also the- 
ories where dark stars accrete matter reaching masses of 
up to 10 7 Mq, w hich would improv e prospects for observa- 
tion significant ly (|Freese et al . 2010; Zackrisson et al. 2010b: 
line et alJl20T3 ). 



5 SUMMARY 

In this paper, we have investigated the prospects of detect- 
ing high-redshift population III stars using the JWST. While 
isolated population III stars or small clusters of population 
III stars at z > 2 will be intrinsically too faint for detection, 
a foreground galaxy cluster acting as a gravitational lens 
could lift them above the detection threshold. We found that 
the detectability of single stars depends strongly on stellar 
mass, with more massive stars more easily detectable. With 
gravitational lensing we find for a 60 Mq star at z ~ 7.5 a 
minimum Logio SFR of -3.4, considering a very high Lya 
escape fraction of 0.5. For a 300 Mq star we find a corre- 
sponding minimum Logio SFR of -4.4. Since our comparison 
simulation model at that redshift indicates a Logio SFR of 
-4.3 this would actually imply a detection of ~ 1 population 
III star. So it turns out that unless the typical mass of pop- 
ulation III stars approaches 300 Mq, their surface number 
density will be too low for detection. 
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